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We reexamine the problem of delocalization of two-dimensional electrons in the presence of random 
magnetic field. By introducing spatial correlations among random fluxes, a well-defined metal- 
insulator transition characterized by a two-branch scaling of conductance has been demonstrated 
numerically. Critical conductance is found non-universal with a value around e^/h. Interest- 
ing connections of this system with the recently observed _B = two-dimensional metallic phase 
(Kravchenko et at, Phys. Rev. B 50, 8039 (1994)) are also discussed. 
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Whether two-dimensional (2D) electrons can become 
delocalized in the presence of random magnetic field 
(RMF) is still controversial. This is a very important 
issue related to many interesting systems, like half-filled 
quantum Hall effect (QHE) [0J^, gauge-field descrip- 
tion of high- Tc superconductor and so on. By us- 
ing the standard transfer-matrix method Q, a number 
of numerical calculations (^-|^] have been performed for a 
non-interacting 2D electron system subject to spatially- 
uncorrelatedKMF . The results indicate that electrons are 
always localized near the band edge, while there is a dra- 
matic enhancement of localization length as one moves 
towards the band center. However, the interpretation 
of the latter is rather confiicting, ranging from that all 
states are still localized with an extremely large lo- 
calization length close to the band center to the exis- 
tence of a critical region Q with divergent localization 
length. Even if a critical region characterized by wave- 
functions with fractional dimensionality could exist 
here, a metallic phase seems being ruled out by those nu- 
merical calculations since a two-branch scaling as a hall- 
mark for metal-insulator transition (MIT) has never been 
found. Analytically, while the study based on a pertur- 
bative nonlinear sigma model approach pointed |^ to the 
localization of all states, the existence of extended states 
was shown possible in the presence of a long-range 
logarithmic interaction of the topological density (due to 
fiuctuating Hall conductance |11]), which is supported by 
direct numerical calculations |12,|l^ of topological Chern 
number for the case of spatially-uncorrelated RMF with 
reduced field strength. 

In contrast to spatially-uncorrelated RMF, however, 
magnetic fiux fiuctuations in realistic systems [|l]-^ may 
be much more smooth with finite-range spatial correla- 
tions. Such a smoothness can significantly reduce the 
random scattering effects while still retain the delocal- 
ization effect |p^-|l3| introduced by magnetic fluxes. In 
this paper, we demonstrate numerically for the first 
time the existence of MIT which is characterized by a 
two-branch scaling of conductance in the presence of 
spatially-correlated RMF. The critical conductance itself 
is non-universal, with its value around e^//i which gen- 



erally increases as the Fermi energy shifts towards the 
band center. With much reduced error bar, the present 
numerical algorithm is also applied to an uncorrelated 
(white noise limit) RMF case and the results unambigu- 
ously show that all states are localized without a critical 
region at strong strength of RMF. Possible connections 
of the present RMF system to the zero-magnetic-field 
(B=0) 2D metal [Q are also discussed at the end of the 
paper. 

We consider a tight-binding lattice model of noninter- 
acting electrons under RMF. The Hamiltonian is defined 
as follows: 
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Here cj is a fermionic creation operator, and < ij > 
refers to two nearest neighboring sites. Wi is an uncorre- 
lated random potential (white noise limit) with strength 
\wi\ < W. A magnetic fiux per plaquette is given as 
4'W — X]n where the summation runs over four links 
around a plaquette labeled by k. We are interested in the 
case where (j){k) at different k's is correlated which can 
be generated in the following way: 
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where Rk (Ri) denotes the spatial position of a given 
plaquette k{i). ho and A/ are the characteristic strength 
and correlation length scale of RMF, respectively, fi is a 
random number distributing uniformly between (-!,+!). 

We employ the following numerical algorithm to cal- 
culate the longitudinal conductance Gxx- Based on the 
Landauer formula, Gxx for a square sample Af ^ L x L 
can be determined as a summation over contributions 
from all the Lyapunov exponents of the Hermitian trans- 
fer matrix product T+T HJl^. To reduce the bound- 
ary effect of a finite-size system, we connect M different 
square samples together to form a very long stripe along 
X direction [of size L x (LM)]. Typically M is chosen 
to be larger than 5000 even for the largest sample size 
(L — 200) in this work. In this way, the statistical error 
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bar is significantly reduced in our results (about 1.5%). 
In most of earlier numerical calculations, finite-size local- 
ization length has been computed where the statistical 
fluctuation is usually quite big (especially near the band 
center) as compared to a direct calculation of the finite- 
size longitudinal conductance in the present algorithm. 

As a test, we have first re-studied the case in which 
the flux (/)(A:) is randomly distributed between —it to tt 
without spatial correlations - the situation investigated 
previously as mentioned at the beginning of the pa- 
per. We find that Gxx monotonically decreases with the 
sample size L at all strengths of the on-site disorders: 
from W = to W = A, and is extrapolated to zero at 
large sample-size limit as shown in Fig. 1 at a fixed Fermi 
energy Ef = —1. In the insert of Fig. 1, Gxx is shown as 
a function of the disorder strength W at different sam- 
ple sizes: L = 24, 80, and 200, which shows that even 
at = the conductance monotonically decreases with 
the increase of L, indicating that the dominant role of 
the random flux here is similar to the random potential 
in causing localization of electrons. The one-parameter 
scaling of Gxx can be obtained by choosing a scaling vari- 
able ^ at each random potential W. As plotted in Fig. 
2, all data can be then collapsed onto a single curve of 
L/^, in which ^ is given in the insert of Fig. 2. Clearly 
^ is always finite although it becomes extremely large at 
weak disorder limit. This is consistent with the conclu- 
sion [^,0 that electrons are all localized and excludes the 
possibility of a critical region as the error bar in our 
calculation is much less than the variation of the con- 
ductance itself. Notice that in weak-disorder limit ^ may 
no longer be interpreted as localization length which 
characterizes an exponential decay of conductance with 
sample size at strong localized region. 

Now let us focus on RMF with smooth spatial corre- 
lations as defined in (2). With the correlation length 
Xf = 5.0 (the lattice constant as the unit) and flux 
strength ho = 1, Gxx as a function of disorder strength W 
is computed at a given Fermi energy Ef = —1 as shown in 
Fig. 3. Curves at different sample sizes (L = 16 - 200) all 
cross at a fixed-point W — Wc, which is independent of 
lattice size L within the statistical error bars. It is qual- 
itatively different from the behavior of Gxx in spatially- 
uncorrelated RMF case discussed above. At W > Wc, 
Gxx continuously decreases with the increase of the sam- 
ple size, which can be extrapolated to zero at large L 
limit, corresponding to insulating phase. On the other 
hand, at 14^ < Wc, Gxx monotonically increases with lat- 
tice sizes like a typical metallic behavior. The insert of 
Fig. 3 shows the critical conductance Gc (corresponding 
to W — Wc) at different Fermi energies and /iq's. The 
data of Gxx in Fig- 3 can be collapsed onto a two-branch 
curve as a function of scaling variable L/^ as shown in 
Fig. 4 for > Wc and W < Wc, respectively. The 
insert of Fig. 4 shows the scaling variable ^ vs. W which 
diverges at the critical point Wc- In the metallic phase 



at < Wc, Gxx can be approximately fitted by the fol- 
lowing form: Gxx = Gs — co*exp{—L/^o). Here Gg is the 
saturated conductance at L ^ cx), which is non-universal 
and depends on the disorder strength W as well as the 
correlation length A/ of random fluxes. 

The introduction of spatial correlations in random 
fluxes is crucial for such a metal-insulator transition. We 
also found a well-defined MIT at an even shorter cor- 
relation length: A/ — 2.0. But the larger A/ is, the 
stronger the metallic behavior becomes with a larger 
saturated conductance. The previously discussed RMF 
in white noise limit may only belong to a very special 
case in which the localization effect of strong random- 
ness of fluxes overwrites the delocalization effect of the 
same fluxes. We would like to point out that even in 
such an uncorrelated random flux case, the delocaliza- 
tion may be still enhanced if one reduces the strength of 
RMF. Earlier topological Chern number calculations |l^ ] 
clearly indicates a delocalization transition as the max- 
imum strength of (/)(/c) is reduced to around tt/2. We 
have computed the conductance in this case using the 
present method at much larger sample sizes and indeed 
found a slight increase of the conductance with sample 
size at W < Wc, which is opposite to strong random 
flux limit where conductance always decreases with the 
increase of sample size (Fig. 1), although a two-branch 
scaling curves here is not as clear-cut as in the spatially- 
correlated RMF case shown in Fig. 4. 

As mentioned above, the critical conductance Gc varies 
from Q.be^/h to around 2e^/ft, as the Fermi energy shifts 
from the band edge towards band center (the insert of 
Fig. 3). It is interesting to note that Gc obtained here 
is in the same range as the experimental data found in 
recent B = 2D MIT system Q. In the following, 
we would like to point out a possible deeper connection 
between the two systems. 

In a recent experiment ||l^ in p-type GaAs/AlGaAs 
heterostructure, the evolution of delocalized states was 
studied continuously from the QHE regime at strong 
magnetic field to zero field limit where the B — MIT 
is recovered. The authors found that the critical den- 
sity of the lowest extended level in QHE regime flattens 
out, instead of floating up towards infinity, as magnetic 
field is reduced and can be extrapolated to the critical 
density of i? = MIT in such a material. Similar result 
has been also observed in Si-MOSFET samples ||l^,|l8|. 
At first sight, it is tempting to think that the lowest ex- 
tended level of QHE somehow survives at £? = 0, but 
physically it does not make much sense because QHE ex- 
tended states carry quantized Hall conductance known 
as Chern number whereas at i? = the total Hall con- 
ductance must be zero without time-reversal symmetry- 
breaking. In fact, experiments indicated JT^ that be- 
fore B vanishes, extended levels of the QHE may al- 
ready merge with a different kind of extended level (called 
QHE/Insulator boundary in Ref, Q) which carries an 
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opposite sign of Hall conductance. Theoretically, it has 
been previously found ||l^ that QHE extended states in- 
deed can be mixed with some boundary extended level 
moving down from high-energy side at strong disorder or 
weak magnetic field limit which carries negative Chern 
number in a lattice model. When those extended states 
with different signs of Chern numbers mix together at 
weak magnetic field limit, there could be two conse- 
quences: one is that no states will eventually carry non- 
zero Chern number due to the cancellation such that all 
of them become localized. This is what happens in non- 
interacting system The second possibility is that 
individual states may still carry nonzero Chern numbers 
and form a delocalized region even though the average 
Hall conductance still vanishes at i? = 0. Such a system 
is then physically related to the RMF system where the 
delocalization mechanism is also due to the fluctuating 
Hall conductance [|T^-|T3|. Below we give a heuristic ar- 
gument how a strong Coulomb interaction may lead to 
such a realization. 

At strong Coulomb interaction with ^ f (here rg 
is the ratio of the strength of the Coulomb interaction 
over the Fermi energy jlj]), the 2D electron state is very 
close to a Wigner glass phase where the low-lying spin 
degrees of freedom may be described by an effective spin 
Hamiltonian Hg given in Ref. ||2^. The low- lying charge 
degrees of freedom may be regarded as "defects" which 
can hop on the "lattice" governed by a generalized t — J 
like model |20y2l| ]. Based on many studies on the t — J 
model in high- Tc problem, especially the gauge-field de- 
scription 1^, charge carriers moving on a magnetic spin 
background can generally acquire fictitious fluxes. Such 
kind of fluxes usually can be treated as random mag- 
netic fields with some finite-range spatial correlations. 
According to the numerical results presented above, such 
a system indeed can have a MIT at i? = 0. Of course, 
further model study is needed in order to fully explore 
this connection which is beyond the scope of the present 
paper. 

In conclusion, we have numerically demonstrated the 
existence of a metal-insulator transition characterized by 
a two-branch scaling for 2D electrons in the presence of 
spatially- correlated random magnetic fields. In contrast 
to usual three-dimensional metal where the conductance 
scales to infinity, this 2D metal has a saturated non- 
universal conductance. The range of the critical con- 
ductance is very similar to that found in _B = 2D 
metal-insulator transition. We briefly discussed a possi- 
ble connection between a 2D interacting electron system 
at ^ 1 and the spatially-correlated random-magnetic- 
fleld problem based on both experimental and theoretical 
considerations. 
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Fig. 1 The evolution of conductance Gxx (in units of 
e^ /h) with sample width L at differnent disorder strength 
Ws for spatially-uncorrelated RMF case. The insert: 
Gxx as a function of W at different Us. Fermi energy is 
fixed at Ef = —1. 
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Fig. 2. The data of Gxx at different L's and W's all 
collapse onto a scaling curve as a function of L/^. The 
insert: £, versus W. 

Fig. 3 Gxx versus W at different sample sizes {L = 
16(«), 24, 32, 48, 64, 80, 120, 200). Wc is the critical dis- 
order. Fermi energy is chosen at Ef = —1. The insert: 
critical conductance Gc as a function of Fermi energy Ej. 

Fig. 4. Two branch-scaling curve of Gxx as a single 
function of L/^ for different L's and W's. The insert: ^ 
versus W. 
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